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Abstract 

We evaluate the shear viscosity of a pion gas in the relativistic kinetic theory approach. The 
in-medium propagator of the p meson at finite temperature is used to evaluate the tt — tt scattering 
amplitude in the medium. The real and imaginary parts of the self-energy calculated from one-loop 
diagrams are seen to have noticeable effects on the scattering cross-section. The consequences on 
temperature dependence of the shear viscosity evaluated in the Chapman-Enskog and relaxation 
time approximations are studied. 
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I. INTRODUCTION 



One of the most exciting revelations that has emerged from experiments at the Rela- 
tivistic Heavy Ion Collider (RHIC) is the strongly interacting nature jlj of the produced 
matter. It has properties of an almost perfect fluid characterized by a rather small value of 
the shear viscosity over entropy density ratio T]/s, close to the quantum or KSS bound 
This interpretation based on the measured elliptic flow v% of hadrons in terms of viscous 



hydrodynamics however depends sensitively on the value of r]/s. A lot of interest has 
generated leading to quite a few estimates of the transport coefficients of both partonic 
as well as hadronic 



been 
4| 



ll| constituents of strongly interacting matter. In general, the shear 



viscosity coefficient is inversely proportional to the cross-section leading to strongly inter- 
acting systems being less viscous than non-interacting ones, the larger values of the viscous 
coefficients resulting from the possibility of momentum transport over larger distances in 
the latter case. The interaction cross-section thus plays the pivotal role in the relaxation of 
systems towards equilibrium. 

Pions form the most important component in a hadronic system. There have been quite 
a few estimates of the viscosity of a pion gas. Almost all these calculations have been 
performed in the relativistic kinetic theory approach in which a crucial input is the irir 
differential cross-section. In jf], 9], the scattering am plit ude was estimated from the lowest 



order Lagrangian of chiral perturbation theory [12|, [13( and in |7| a unitarized amplitude 
was employed for a meson gas. Phenomenological amplitudes obtained from fits to phase 
shift data have been employed in [h3| in view of the fact that the tttt cross-section estimated 
from lowest order chiral perturbation theory is known to deviate from the experimental data 
beyond centre of mass energy of 500 MeV primarily due to the p pole which dominates the 
cross-section in the energy region between 500-1000 MeV. 

The spectral modification of the p in hot and dense matter is known to play the most 



important role in explaining the invariant mass spectra o 



ion collisions in the region below the nominal p mass 14 



lepton pairs from relativistic heavy 



161 ] . The lepton pair production in 



this region is related by vector meson dominance to two-pion annihilation proceeding through 
p-exchange. Moreover, medium modification of the p spectral function incorporated in the 



evaluation of elliptic flow of lepton pairs of low invariant mass have been recently shown 



17| 



to result in values close to those of hadrons of similar mass, ascertaining the role of Vi of 
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hadrons in describing the early stages of heavy ion collisions. It is natural to extend this 
scenario to the evaluation of the transport coefficients, the temperature dependence of which 
is particularly relevant for the electromagnetic probes which are emitted at all stages of the 
collision spanning a range temperatures. Investigations with temperature dependent titt 
cross-section involving a exchange have recently been addressed in [l8|. Our intention in 
this work is to observe the effect of the in-medium spectral function of the p on the shear 
viscosity of a pion gas. 

The first task is to obtain the energy dependent irir cross-section using a phenomenological 
approach which is close to the experimental value and at the same time is theoretically 
amenable to the incorporation of medium effects. To this end we consider the scattering 
to proceed via p exchange for which the invariant amplitude is evaluated using effective 
interactions. Medium effects on p propagation are introduced through one-loop self-energies 
to obtain the modified cross-section at finite temperature. These are discussed in section 
II. The corresponding effect on the shear viscosity of a pion gas is the subject of section 
III where it is evaluated using the non-relativistic, Chapman- Enskog and relaxation time 
approximations. We end with a summary in section IV. Details of the calculation of the p 
self-energy in the medium has been summarized in the Appendix. 



II. THE vrvr CROSS-SECTION WITH MEDIUM EFFECTS 

We evaluate the invariant amplitude for -kit scattering using an effective Lagrangian in 
which the coupling of the p meson to the pions is introduced through the gauge covariant 
derivative of the pion field operator to obtain [3] 

C^= l -^Tr[V^[d^M\ (1) 

where Tr indicates trace in SU(2) space. The matrix $ collects the pion fields in the form 

vr v / 2vr H 



x/2tt- - !l 



and V M collects the p meson fields analogously. The coupling constant 

<?p=6.05 is fixed from the p — > tttt decay width. This interaction leads to tttc scattering 
diagrams with p exchange in the s, t and u channels. In these calculations we modify the 
p propagator Dffi = {—g^ u + QfiQu / 'ntf,) I '(q 2 — m<2 p + ie) replacing ie with im p T p (s) where the 

2 

two-pion decay width is T p (s) = j$^(s — 4m^) 3//2 . This is done only for s-channel p-exchange 
diagrams which contribute only in the case of total isospin 7 = 1. 
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FIG. 1: The irir cross-section as a function of centre of mass energy. The dashed line indicates the 
cross-section obtained using eq. [2] which agrees well with the experimental values shown by filled 
circles. The dash-dotted and solid lines depict the in- medium cross-section for mr and multi-pion 
loops respectively in the p self-energy evaluated at T=160 MeV. 



In order to describe 7nr scattering at low energies it is essential to also include a-exchange 
diagrams. We use the interaction C a7T7T = \g a m a Ti ■ na with g a = 2.5 to calculate the 
scattering amplitudes where as before we have introduced the a width in the s-channel 
diagram which now appear only for 1 = 0. The values m a = 450 MeV and IV = 550 MeV 



that we use are in conformity with estimates in [20]. Going over from the charge to the 
isospin basis it is easy to write down the matrix elements for all the possible values of the 
total isospin of the two pions. They are, 



M I=0 
Mi=i 
M I=2 



St 



s — u s 

-z + — 



t 



t 



rn 



u — mz 



p - ■■~p 
2(t - u) 
s - m 2 p + im p V p (s) 

u — s t — s 
k + 



, 2 2 
+ 9a m a 



1 



1 



s — ml + im a T a t — ml 



+ 



t - s 



u — s 



u — m 



u — mt 



i 2 2 
+ 9a m a 



t 



m 



2 2 
9a m a 



u — m" D 
1 



u — mt 



t — m 2 „ 



u 



mt 



(2) 



The differential cross-section is then obtained from ^ = \M.\ 2 / QA-k 2 s where the isospin 



averaged amplitude is given by \A4\ 2 = | 5^(21 + 1)|A^/| 2 . Ignoring the 1 = 2 contribution, 
the integrated cross-section (with an additional factor of 1/2 for identical particles) plotted 
as a function of the centre of mass energy in Fig. [1] (dashed line) is seen to agree reasonably 
well with the one taken from Ref. 2jJ (shown by filled circles) which corresponds to a 
resonance saturation parametrization of isoscalar and isovector phase shifts obtained from 
various empirical data involving the tctt system. We will use this theoretical estimate as a 



benchmark to study medium effects. 

The effect of the medium on p propagation is quantified through its self-energy. The 
standard procedure is to evaluate this quantity by perturbative methods using effective 
interactions and then obtain the exact propagator using the Dyson equation, depicted pic- 



torially in Fig. [2J We now reco 
have been explicitly given in 22 



lect the main steps in this scheme omitting details which 
time formulation of thermal field theory, all 



23| . In the rea 



24( which can be diagonalized. The diagonal 



two-point functions assume a 2 x 2 matrix form 
components also obey the Dyson equation 25j by means of which the full propagator D^ v 
is obtained as 

D,u = D^+D^W"D Xu (3) 



where D^J is the vacuum propagator for the p meson and n CTA is the self energy function 



obtained from one-loop diagrams shown in Fig. [2j Following 
self-energy in terms of longitudinal and transverse parts 



22. 



24| we write the in-medium 



(4) 



where P^ v and Q v ,, are the transverse and longitudinal projection tensors respectively. These 
are defined as 22] 

,2 



P, 



fj,v 



9fJ.V T 9 _ 

q z q 



^u^Uy, Un = Up - (u ■ q)q ll /q 



and 



Q - <^ - - - - - 



z^-u^Uu, q = (u- q) -q 
q z 



where is four velocity of the thermal bath. It is easy to see that 



Pnv + Q^u/q 2 = -g» v + q^qu/q 2 



(5) 



(6) 



(?) 



Note that while P and Q are four-dimensionally transverse, P is also three-dimensionally 
transverse while Q is longitudinal. Solving (J3J), the exact p propagator is obtained as 



Dnu{qo,q) 



P 



[IV 



Qpiv/q 1 



+ 



q 2 — m 2 — IT 11 q 2 — m 2 — q 2 Yi L m 2 q 2 
where I1 T and U L can be obtained from the relations 



(8) 



it =- (n^ + ^n 00 ), n 

q l 



q z 



(9) 
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FIG. 2: The exact /> propagator with ir — h loop diagrams for h = ir,u),hx,ai mesons. 

As shown in [22j, the three-momentum dependence of the p self-energy is not substantial 
for the present case and we can replace n T and q 2 U L in the above expression by a self-energy 
function which is averaged over polarization. Defining this as 

n = ^(2n T + g 2 n L ) (10) 

and neglecting the non-pole piece in (JHJ), the in-medium propagator can be written as 

— g ^ _|_ q Qi//Q^ 

D^(qo,q) = — 2 2 — tJTF? — V, -t tt^ ^ ' ( U ) 

q 2 - m l p - Rell(go, q) + umlI(go, q) 

The real part of the self-energy modifies the pole position and the imaginary part embodies 
the effect of collisions and decay processes by means of which the p is lost or gained in the 
medium. Using interactions from chiral perturbation theory, the one-loop self energy was 



recently calculated 22]. Some details of the calculation are provided in the Appendix. The 
imaginary parts for tttt, ttu, it hi and na\ loops were obtained from the discontinuities of 
the self-energy in the complex energy plane. While for the tttt loop, the contribution at the 
nominal p pole comes from the unitary cut, the Landau type discontinuity is responsible for 
contributions from loops with heavier particles. The mesons ou, hi and a\ all have negative G- 



parity and have substantial 3ir and pn decay widths 26( . The self-energies containing these 



unstable particles in the loop graphs have thus been folded with their spectral functions as 



shown in 23] . The contributions from the loops with heavy mesons may then be considered 
as a multi-pion contribution to the p self-energy. The cross-section obtained by using the 
in-medium p-propagator (fTTj) in place of the vacuum propagator Dffj in the evaluation of the 
amplitudes is shown in Fig. [TJ We observe a small suppression of the peak for the tttt loop 
and a larger effect when all the loops (indicated by multi-pion) are considered accompanied 
by a small shift in its position. This is due to the temperature dependence of the real and 
imaginary parts of the self-energy and is manifested as the modified spectral function of 
the p meson. Extension to the case of finite baryon density can be made using the results 
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FIG. 3: The classical estimate of the shear viscosity of a pion gas as a function of temperature. 
The inset shows the behaviour at low temperature. The dashed line with circles depicts the s/T 
behaviour at low T. The other curves in the inset are the continuations of those shown in the main 
figure to lower temperatures. 



in 23J]. Similar reduction of the tttt cross-section for a hot and dense system were also 



obtained earlier in 127] . 



III. THE SHEAR VISCOSITY OF A PION GAS 



A first insight into the medium effects on the shear viscosity may be obtained by em- 
ploying the classical estimate of the latter for a dilute gas of pions in terms of the mean free 
path A = 1/ na and is given by 



(12) 



where the thermally averaged pion momentum and cross-section for the process 7r(p) + 
ir(k) — > 7r(p') + ir(k') are given respectively by 

Jd 3 p \p\f(E p ) 



p(T) 



fd*pf(E p 



_ = Jd 3 p d 3 k f(E p )f(E k )a(s)(l + f(E p/ ))(l + f(E k ,)) 
° ~ fdZpd*kf(E p )f(E k )(l + f(E p ,))(l + f(E k ,)) 



(13) 



where f(E) denotes the Bose-Einstein distribution and E p = ^p 2 + m1 and E k = 
\Jk 2 + m\. Manipulation of the exit channel in the cm. frame leads to E p i = Ey = 
E c .m./2 = y/s/2. We plot the shear viscosity as a function of temperature for a range of 
values relevant for the hadronic phase of heavy ion collisions in Fig. [3j The behaviour of 



7] in this picture is determined by the interplay of the temperature dependence of p (which 
essentially behaves as yT since p 2 /2m ~ kT ) and that of a. A substantial effect of the 
reduced in-medium cross-section can be seen in this approximate estimate of 77. In spite of 
possible differences in the averaging procedure, the values of r] for the vacuum case lie within 
similar range with those in [lOj. The behaviour at low temperatures is shown separately in 
the inset of Fig. [3j Since the energy of scattering decreases with T as shown in [hJ, the 
cross-section tends to a constant as the temperature is lowered (Fig. [1]). The shear viscosity 
is then proportional to y/T, a situation depicted by the dashed line with circles (for arbitrary 
proportionality constant). This \/T behaviour of the shear viscosity at low energies has been 
shown analytically using a hard sphere tttt cross-section in [6|. We also note in the inset that 
all the curves actually converge as the temperature is lowered showing that medium effects 
reduce with pion density. 

In quantum field theory, the shear viscosity is defined using linear response in terms of 
the correlation function of the (spatial component) of the energy-momentum tensor which 



can be obtained from the Lagrangian density of the system |28|, [29|. But as argued in |3CH. 13 1 1 
a perturbative (diagrammatic) evaluation involves a sum over an infinite set of complicated 
diagrams in order to obtain the leading behaviour even in weakly coupled systems. A more 
efficient approach is to obtain these by solving the linearized transport equation. 

The relativistic transport equation for the phase space distribution f(x,p) of a pion gas 



is written as 



32| 



p^f(x,p) = C[f] (14) 

where C[f] is known as the collision term. For binary elastic collisions p + k — )■ p' + k! the 
Uehling-Uhlenbeck collision term is given by 

C[f] = JdT k dT p , dr k ,[f(x,p')f{x,k'){l + f(x,p)}{l + f(x,k)} 

-f(x,p)f(x,k){l + f(x,p')}{l + f(x,k')}} W (15) 

where dT a = lr f,i and the collision rate W is defined as 

W= S -^(2n)^(p+k-p>-k') 

which includes an explicit factor of 1/2 for indistinguishable particles. 

A system in local equilibrium is described by ideal fluid dynamics and departures from 
this state can be dealt with using dissipative dynamics. The correspondence between non- 
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equilibrium kinetic theory and viscous hydrodynamics can be studied by considering small 



departures from equilibrium. Following 



33|, in the first Chapman-Enskog approximation 



one takes the non-equilibrium distribution as 

f{x,p) = fW(x,p) + f {0 \x,p)[l + f°\x,p)]4>(x,p) (16) 

where f(°>(x,p) is the local equilibrium (Bose) distribution with space-time dependent pa- 
rameters and <p(x,p) is the deviation therefrom. Putting in ( 1T41) the function <p(x,p) is seen 
to satisfy the linearized transport equation 

p^fW(x,p) = -C[(j>] (17) 

where 

£{<p] = /(°) J dv k dv p , dr k ,f(°\x,k){i + ;(°W)}{i + f (0) ^k')} 

[4>{x, p) + (f>(x, k) - <j){x, p) - <p(x, k')] W . (18) 

From general considerations it can be shown that in order to be a solution of this equation, 
the function has to be of the form = Ad^u^ + C P (9V) neglecting terms related to 
thermal conduction. The symbol (• • •) indicates a space-like symmetrized and traceless 
combination. 

To make a connection with the coefficient of viscosity one makes use of the energy- 
momentum tensor given by 

TT = J dT pP ^f(p) . (19) 
Putting in the non-equilibrium distribution function ([TBI) into this equation, T^ v separates 



into ideal and dissipative parts. The latter is then compared with the most general form 34, 
35| for this quantity to arrive at the expression 



7/ 



~ / dT p C, u {p^)f Q \p){\ + /(°>(p)} . (20) 



Obtaining C^ u by solving eq. ([IT]) is the most important component of the calculation. 
Evaluation of the advective derivative of the equilibrium distribution function on the left 
hand side of this equation results in the integral equation 

C[C, U ) = -(p»p„)f i0) (p){l + / (0) (p)}/T (21) 
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FIG. 4: The shear viscosity as a function of temperature in the Chapman-Enskog approximation. 
The upper set of curves were obtained with ipmax ~ 2 in eq. 1241 where the dash-dotted line indicates 
use of the vacuum cross-section and the dashed and solid lines correspond to in-medium cross- 
section for the 7T7T and multi-pion cases respectively. The lower set of curves (with circles) were 
obtained with ipmax = oo. Here, the dash-dotted line (with circles) represents the vacuum case and 
the dashed and solid lines (with circles) correspond to in-medium cross-sections for the tttt and 
multi-pion cases respectively. The difference between the corresponding curves in the upper and 
lower sets is a measure of the uncertainty in r\ due to insufficient information about the cross-section 
for energies more than ~ 1 GeV. 

which can be reduced to algebraic equations by expanding the unknown function in 
terms of orthogonal polynomials. Similar techniques have been used by many authors to 
calculate the transport coefficients [6l4l0l. Il8|. Using e.g. generalized Laguerre polynomials 
(of half-integral order) the first approximation to the shear viscosity obtained by keeping 



the first term in the series is given in the notation of 

10 c 00 

where 



33 



m by 



(22) 



To 
coo 



-10 



S 2 \z) 
h(z)+I 2 (z)+I 3 (z) 



with z = m^/T and 



(23) 



The integrals I a (z) are given by 

N ( z\ r°° r 77 1 dcr t°° 

I*(z) = ° 12 / cosh^sinhV / cosine- — (V>, 6) / d X L a ( X ) 
[b 2 [z)f Jo Jo 2dil Jo 
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where 



2n rw p2z cosh ip cosh x 

d(f> d9 sine . E F — ; G 1V g - M tt (g,6) (24) 

o 7o (e B - l)(e F - l)(e G - l)(e H - 1) 



= z(cosh ip cosh x — sinh -0 sinh x cos 9) 

F = z(cosh ip cosh x — sinh ■?/> sinh x cos 9') 

G = E + 2z sinh ip sinh x cos # 

= F + 2z sinh ip sinh x cos 6 1 ' , (25) 



and the functions N, L, M and 5 are defined as 



Nt(z) = I6z\ N 2 (z) = 16z A , N 3 (z) = 16z 4 /3 , 

L x {x) = sinh 2 x, L 2 (x) = sinh 4 x, L 3 (x) = sinh 6 x , 

M x (0,e) = l-cos 2 9, M 2 (^,6) = cos 2 + cos 2 ^-2cos^cos^cose , 

M 3 (^,e) = [cos 2 #-cos 2 #f and 

oo 

^(^ = E kaK n(kz) . (26) 
fc=l 

Here K n (x) denotes modified Bessel function of order n and cos#' = cos 9 cos — 
sin 9 sin G cos <p ■ 

We plot t] versus T in Fig. H] obtained in the Chapman- Enskog approximation showing 
the effect of the in-medium p propagation in the pion gas. In the upper set of curves the 
upper limit of integration over ip in eq. (124)) is taken to be ~ 2 which corresponds to a center 
of mass energy (E c m ^ = 2m n cosh ip) of ~ 1 GeV for tttt scattering up to which we have 
obtained a fair description of the cross-section using p and a exchange. Here we observe 
~ 10% change at T = 150 MeV due to medium effects compared to the vacuum when all 
the loops in the p self-energy are considered. The effect reduces with temperature to less 
than 5% at 100 MeV and the curves almost merge with each other as T reaches ~ 60 MeV. 
The lower set of curves (with circles) in Fig. H] depict the corresponding set of curves when 
the actual upper limit of ip (which is oo) is used. The difference between the two set of 
curves indicates the uncertainty in the results on account of insufficient information of the 
cross-section at higher energies. 

We also present results for the relaxation time approximation which is the simplest way to 
linearize the transport equation. Here one assumes that f(x,p) goes over to the equilibrium 
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FIG. 5: The shear viscosity as a function of temperature in the relaxation time approximation. 
The dash-dotted and solid lines correspond to the use of in- medium cross-sections in eq. [27] for tttt 
and multi-pion loops respectively. The dashed line represents the vacuum case. 

distribution f^(x,p) as a result of collisions and this takes place over a relaxation time r{p) 
which is the inverse of the collision frequency u(p). The right hand side of eq. (|14p is then 
given by — E p co(p)[f(x,p) — f^°\x,p)]. For collisions of the form n(p) + n(k) — > n(p') + %(k') 
we have 



uj(p) 



dTk 



s(s 



4ml] 



2 .EL 



-f(E k )(i + f(E p ,))(i + f{E k ,)y- / dn^ 



(27) 



23, 



and the shear viscosity is given by 

V = ToT / dV ^ r{p)f{Ep){l + f{Ep)) 



(28) 



We show the temperature dependence of r\ in the relaxation time approximation in Fig. |5j 
The values in this case are lower than that obtained in the Enskog method though the effect 
of the medium is larger. The difference depends largely on the energy dependence of the 



cross-section as noted in 



39|. 



IV. SUMMARY AND OUTLOOK 



To summarize, we have used an effective Lagrangian approach to evaluate the tttt scat- 
tering cross-section in which the temperature dependence enters through the in-medium 
propagator of the p which mediates the interaction between the pions. This has been used 
to evaluate the shear viscosity of a pion gas to the first Chapman- Enskog order. The tem- 
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perature dependence with and without medium effects shows a noticeable difference in this 
approach as well as in the relaxation time and non-relativistic approaches. 

We end by stating that the main focus in this work has been to emphasize the role of 
medium modifications of the cross-section in the evaluation of the transport coefficients. The 
values of r\ arrived at in the present calculation may be modified in a more realistic scenario 



employing model independent approaches in the lines of |40j. In addition, the treatment 
has to be extended to include other mesons and nucleons. Only then can the temperature 
dependent shear viscosity be included in the hydrodynamics instead of constant values used 



411 ] . Such an 



e.g. in the evaluation of electromagnetic spectra from heavy ion collisions 
effort is underway and will be reported in due course. 

V. APPENDIX 



Using standard techniques of thermal field theory, the polarization function of the p 



corresponding to the tt — h{= ir, u, hi, a\) loop diagrams shown in Fig. [2] is given 22|, |23J by 



rJ 3 k 1 



(2tt) 3 4uv^ 



fJ,V 



+ 



g - uJn - w A + ir]e{q Q ) q -u 1T + uj h + ir]e(q ) 

(A.l) 



+ /K)ivr - /K)jyr + -fMNr - a + /k))*t 



q + u 7r -u h + ir]e(q ) q + u w + u h + i7]e(q ) 



where co n = yk 2 + m 2 and Uh = y (q — k) 2 + m\ are the Bose distribution functions. The 
factor contains tensor structures from the vertices and the vector propagators and are 
given by 



2 



N$(q,k) = -4(j|) (B, u + q 2 k 2 A, u ) 



where 



Nfr\q,h) = -2 if)" '(B^-^C^) (A.2) 



A a p{q) = -g a p + q a qp/q 2 , 
B a/ 3(q, k) = q 2 k a k p - q ■ k(q a kf3 + k a qp) + (q ■ k) 2 g a/3 , 

C a p(q, k) = q 4 k a k/3 - q 2 {q ■ k){q a kp + k a q$) + (q ■ k) 2 q a qp. (A. 3) 
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The interaction Lagrangians used to obtain these are given in [22|]. The subscript i(= 1, ..4) 
on in flA.lj) correspond to its values for k = ui n , — ui n , q — Uh, qo + <^h respectively. The 
imaginary part is easily read off from flA.ll) and is given by 

d 3 k 



x 



[iVf {(1 - f(u n ) - f(u h ))5(q u h ) + (f{u n ) - f(uj h ))6(q + u h )} 

+NP{(f(u h ) - fM)6(q u h ) - (1 - fM - f(u h ))8(q + oj h )}](AA) 

where the factors N34 have converted to N± u or N^ u on use of the associated (^-functions. The 
latter actually define the kinematic domains where different scattering and decay processes 
leading to the loss or gain of the p mesons in the medium. The non-vanishing regions 
produce branch cuts in the self-energy function. The first and the fourth terms are non-zero 
for q 2 > {rrih + m^) 2 giving rise to the unitary cut and the second and third terms are 
non-vanishing for q 2 < (nth — m n ) 2 giving rise to the Landau cut. Whereas the unitary cut 
is present in vacuum as well, the Landau cut appears only in the medium. 

Performing the angular integration, the unitary and Landau cut contributions for the 
physically relevant region q , q 2 > are obtained respectively as 

lmU^(q ,q) = -Jjj^j jf dx N» v (x){l + f(u) + f(q - u)} 

for q > \j{m h + m^) 2 + \q\ 2 
Imnr(g ,g) = ^ jyxN^(x){f(oo)-f(q + u) 

for \q\ < q < \J{m h - mj 2 + \q\ 2 (A.5) 

where x is defined in the two cases through u = (R 2 /2q 2 )(q + \q\x) and u = (R 2 /2q 2 )(—q + 
\q\x) with R 2 = q 2 — ra\ + m 2 and v = \Jl — ^S^- The transverse and longitudinal 
components for different graphs can now be obtained using the relations, 

A» = -3, A 00 = % (A.6) 
B> = mW + y , B 00 = -^(1 - x 2 ) (A.7) 

C* = ftrntf - ^), C 00 = ^x 2 . (A.8) 

The real parts are obtained as principal value integrals which remain after removal of the 
imaginary part from eq. flA.l j) . Note that unlike the imaginary part, the real part of the 
self-energy at a given value of go receives contribution from all the terms. 
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To take into account the finite 3tt and pn decay widths (Th) of the unstable mesons in 
the loop, the self-energy is folded with their spectral function, 



u{q > mh) =w h L 



{m h +2T h f 1 

dM 2 -lm 

h J(m h -2r h ) 2 IT 



1 



M 2 — m\ + iMT h (M) 



U(q,M) (A.9) 



with Nh 



(m h +2T h f I 

dM 2 -lm 

(m h -2T h ) 2 71 



M 2 -m 2 +iMT h {M) 
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